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Abstract 

We use the entropy production variational method to associate a one particle distribution function 
to the assumed known energy-momentum and entropy currents describing a relativistic conformal 
fluid. Assuming a simple form for the collision operator we find this one particle distribution 
function explicitly, and show that this method of linking the hydro and kinetic description is a non 
trivial generalization of Grad's ansatz. The resulting constitutive relations are the same as in the 
conformal dissipative type theories discussed in J. Peralta-Ramos and E. Calzetta, Phys. Rev. D 
80, 126002 (2009). Our results may prove useful in the description of freeze-out in ultrarelativistic 
heavy-ion collisions. 
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I. INTRODUCTION 



The relationship between the kinetic and hydro dynamic descriptions of a dissipative rel- 
ativistic system is a long standing puzzle, because while the relativistic Boltzmann equation 



is well understood 11. 2|, its obvious match, namely the Eckart and Landau-Lifshitz theories 

n 

[3[, are plagued with causality and stability problems. Solving this puzzle has acquired a 
certain urgency, because being able to associate a one particle distribution function to known 
hydrodynamic currents is an essential step in describing freeze-out in hydrodynamic models 
of heavy ion collisions (see for instance 4tS|)- 

There is a long history in the development of different formalisms to derive hydrodynamics 
from kinetic theory, the most well-known methods being Grad's method of moments and the 
Chapmann- Enskog expansion [l|, 2), [^15]. For some recent theoretical developments see 



16]. 



Here, we will use the entropy production variational methods (EPVM) to provide a linkage 
between the kinetic and hydrodynamic description; for a review see jl^]. The idea of entropy 
production variational methods is to find the particle distribution function which extremizes 
the entropy production, providing in the process a means of closing the infinite chain of 
hydrodynamic equations (i.e. providing a closure) 171. Il8j]. 

In order to carry out this derivation explicitly we assume, for simplicity, a conformally in- 
variant theory and use a linear collision operator. This operator satisfies energy-momentum 
conservation and guarantees the H theorem, ensuring that the resulting theory exactly satis- 
fies the Second Law. The size of the collision term is determined by a dimensionful parameter 
r with the physical meaning of a relaxation time. There is a natural expansion of the solution 
in powers of r. We shall find explicitly the one-particle distribution function singled out by 
the EPVM to second order in r. 

The most important result from this paper is the comparison of the distribution function 
picked up by the EPVM and the better known solution provided by Grad's ansatz. We shall 
see that they differ in two main ways. First, the EPVM is not tied up to a gradient expansion 
in the hydrodynamical variables, and so it becomes an attractive choice in situations where 
those gradients are expected to be large. Second, the Grad ansatz is usually applied in 
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conjunction with simplifying assumptions regarding the Boltzmann equation, to the effect 
that the resulting correction to the distribution function is expressed solely in terms of the 
shear viscosity and the equilibrium energy density and pressure; it bears no memory of the 
collision term, unless through a global factor. The correction to the distribution function 
derived from the EPVM, on the other hand, depends in an essential way on the collision 
term. Indeed this result could be seen as a way to associate a simple kinetic equation to 
nontrivial freeze-out prescriptions derived from more fundamental physics and discussed in 
Q (see also fl). 

The EPVM provides a prescription to associate a one-particle distribution function to 
given macroscopic currents, but gives no information on those currents or their further 
time evolution. To fill this gap, we show that the EPVM may be cast in the framework of 
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conformal divergence-type theorys (DTTs). These theories were discussed in [19|, |20j|. In this 
way we determine the dynamics of the energy-momentum and entropy currents; the former is 
conserved, while entropy production is equated to its extremal value. The association with 
DTTs, moreover, affords a simple way of investigating the causality and thermodynamic 
consistency of the hydrodynamics associated to the EPVM. 

As a byproduct, the derivation given here provides a novel kinetic interpretation of the 
nonequilibrium tensor of DTTs, as a Lagrange multiplier enforcing the energy-momentum 
constraint when extremizing the entropy production. 

Divergence-type theories 2l| are hydrodynamic theories which are based on extending 
the set of hydrodynamic variables used to describe a nonequilibrium system to include a 
traceless tensor which vanish in equilibrium. These are exact hydrodynamic theories in the 
sense that they are not based on gradient expansions, and therefore can describe situations 
with large gradients (i.e. shock-waves) in which the so-called second-order theories (SOTs) 
0, [22, 23 1 are known to fail. In 19] we have developed a quadratic DTT for a conformal 
field, which reproduces the SOT developed in 23J when velocity gradients are small. We 
then applied it in [20| to describe the evolution of the fireball created in ultrarelativistic 
heavy-ion collisions. 

The paper is organized as follows. In section [XT] we develop the EPVM approach. We 
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obtain the variational equation for the entropy production, and then we describe the linear 
collision term used here and solve the variational equation perturbately up to second order 
in the relaxation time. In section ITTT1 we compare those results to Grad's quadratic ansatz. 
In section [IV] we overview DTTs and give a brief summary of the main results for confor- 
mal fields, and make explicit the connection between the EPVM and DTTs. Finally, we 
summarize our results in section IVl 



II. ENTROPY PRODUCTION: VARIATIONAL METHOD 

In this section we set up the problem within the context of kinetic theory and derive the 
variational equation to be solved perturbately. We then specify the collision operator and 
solve the variational problem to second order in the relaxation time. 



A. The variational equations 

We consider a relativistic kinetic theory in flat space-time I2I, \d- We use signature 
(—, + + +). The distribution function f(x^,p^) determines the energy momentum tensor 
T^ u and the entropy current 



T^ = J Dpp^f (1) 

and 

S» = J Dpp»{(l + f)\n[l + f]-f\n[f]} (2) 

with 

Dp = -£*e {f) P (V) (3) 

(2tt) 

where p is the density of states and — p 2 = (p ) 2 — p 2 . We assume a conformal theory where 
p{—p 2 ) oc 5 (— p 2 ). The distribution function obeys a Boltzmann-like equation, which we 
write in compact form as 



P»U = hoi [f] ■ 
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(4) 



Energy- momentum conservation implies the constraint 



j Dpp^I col [f] = (5) 
identically in /i and in /, while entropy creation reads 

S^ = J ^/^[/Jlnfl + r 1 ] . (6) 
Given a vector (3^ we define the thermal distribution 



with (3p = It is convenient to parameterize deviations from thermal equilibrium as 

follows 



/ = / [l + (l + / ) x ] . (8) 

Then 

where Tq V corresponds to a perfect fluid, and 

IF- = J Dppp^x (10) 

with 

D f) p = Dpf (l + f ) . (11) 
It is convenient to introduce the notation 



(•••)= / D p p{-..). (12) 



In this notation IP" = (p^p u x) 
We have 



S" = - (3 U W V + S$ = - p u T^ v + S£ (13) 
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where p — p/3 is the pressure and 

Sf = J ^^{(/ - 1 + X )ln[l + /oX]-((l + /o)" 1 + x)ln[l + (l + /o)x]} ■ (14) 



Similarly 



S^ = J D P I col [/]|ln [l + / - 1 ]+ In 



l + /oX 



l + (l + /o)x 



(15) 



But the first term integrates identically to zero, so 



D pP b [x] In 



1 + ZoX 



l + (l + /o)x 



(16) 



where 



Mx] 



hoi [/o[l + (l + /o)x]] 



(17) 



h (i + k) 

By the way, we notice two identities. By taking a variation of the energy-momentum 
conservation constraint, we get 



Dpqq 



u s Ji M (?) 



(18) 



On the other hand, observe that an infinitesimal 5x = 8(3^ is just a shift in (3, and therefore 
1/3 [Sx] ip) = identically in p. Expanding to first order in 5(3 we get 



8Jl M (P) 

^x [q\ 



q 



. 



(19) 



In other words, p^ is both a right and left null eigenvector of 51 p [x] (q) / 5x [p] ■ 

Below we shall restrict ourselves to the case of linear collision terms, these being essentially 
the only ones for which a closed-form solution may be found. The form of the entropy 
production suggests that to conform to the H theorem, such a functional must be linear not 
in x but in the new variable 



C = ln 

This relation may be inverted to yield 



l + (l + /o)x 
1 + ZoX 



(20) 
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1 + (1 + / )X oC 



l + /oX 

so 



(21) 



X - 1 - Jo 1) (22) 
or else, expanding to second order in ( 

X = C + ^(l + 2/ )C 2 • (23) 

Suppose now we wish to find the distribution function that extremizes S 1 ^ given the values 
T^ u and of the energy momentum tensor and entropy current. Choosing a suitable and 
decomposing both the kinetic theory and the hydrodynamic currents as above, we end up 
solving the variational problem 

^ K - A„S? - A^] = (24) 

where A^ and X^ u are Lagrange multipliers enforcing the constraints. Observe that A Ml/ is 
dimensionless while A M has dimensions of temperature. Below we shall restrict ourselves to 
the case A M = 0. 



B. The collision integral 

It is clear that to actually solve for ( we need to know something about the collision 



operator. In this section we shall investigate the structure of linear operators 

ay 

In principle one would like to write Ip [(] = —F£(p) /r, where r is the relaxation time. 
Observe that F has dimensions of temperature, and r has dimensions of time, namely inverse 
temperature. To avoid picking up a preferred direction in the rest frame, it is natural to 
request that F = F [co p ], where u p = —u^. However, these restrictions are not sufficient, 
since such a kinetic equation would violate the energy-momentum conservation constraints. 
To preserve the momentum constraints we introduce a projection operator Q such that 

(p»Q [/]> = (25) 
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for any /, but Q[f] — f if (p^f) = 0. The notation () is defined in eq. ffT2~j) . In the rest 
frame, we write 



Q\f] = f- j-o\ H K/> + W (*/>] (26) 

where we exploit the fact that for a conformal theory (ppi) = 8 V ' (u>T) /3. Now we write the 
collision integral as 

h be] = [FQ [C]] (27) 
Suppose we wish to solve the equation 

Q [FQ [g]] = f (28) 
where (p^f) = 0. Then, in the rest frame 

FQ [g] = f- Au p - Brf . (29) 
The constants A and Bi must enforce the integrability conditions 

(^(f-Aco p - Brf) ) = (/ - Au; p - V) ) = • (30) 

Therefore 

i^p/F) 

and 

<? = 1 (/ - Au p - 5^) . (33) 
We will make use of this properties in what follows. 
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The linear collision operator used here is quite general, and it is interesting to note that 



for F = T we recover Marie's relativistic generalization of the BGK model [24|, while for 
F = u p we get the Anderson- Witting model 25j (see 26| for a comparison of these kinetic 
models to Israel- Stewart formalism). 



C. Perturbative solution 

The structure of the collision term suggests we seek a solution as an expansion in powers 
of the relaxation time r. We shall consider the solution up to second order. 
The equation we wish to solve is 



Q [FQ [(}] = rX^tf [! + (! + 2/ ) C] ■ 



We expand 



and 



Then we find the equations 



C = Ci + C 2 



\ = \(o) + 



(34) 



(35) 
(36) 



Q [FQ [Ci]] = rA^Vp" 

and 

Q \FQ M = rff [A« + (1 + 2/ ) A^Ci] 
and the integrability conditions 



(37) 
(38) 



and 



A(°J ((1 + 2/ ) Ci (P)PW> + Ag) (pW> = 
The correction to the one particle distibution function reads 



(39) 
(40) 



X = Ci + C2 + ^(l + 2/o)C? 
Therefore for the energy momentum tensor we shall find 



(41) 



where 



and 



The entropy flux 



nf = / % P y 



C 2 + ^(l + 2/ )Ci 2 



(42) 
(43) 
(44) 



(45) 



so there is no first order correction. We shall only consider the first nonvanishing contribution 
to the entropy creation 



SI. = A<°> 



(46) 



which is already quadratic in deviations from equilibrium. 
We will now go over to calculate Ci and £2- 



1. First order solution 

Let us now consider the first order solution for the specific collision term introduced above. 
To make things simpler, we shall go to the rest frame, where the integrability conditions 
become 

Ao°o + = (47) 

and 

Ag> = . (48) 
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For simplicity and without loss of generality we shall assume that Aq°j = A- '' 1 = A^ = 0. 
The solution reads 



T 



Ci = j^W ■ (49) 

In our particular case, A = B { = 0. 

If we use this to compute we get IT? = IT? 1 = 0. To compute IT/ recall that for any 
G[u>] 

(G [uj p ] pV'pV) = ^fA^A \s^5 kl + 5 lk 5 lj + 5 a 5 jk ] (50) 

15 

so 

^ = %(f)* m ( 51 ) 

which is traceless. 

The lowest order nontrivial contribution to the entropy flux is 

s? = i?(^H° ,A< °"' (52) 

and S\ — in the rest frame. Similarly 

Si = I (f) ^ m ■ (53) 



2. Second order solution 

We now use the first order solution to investigate the second order one. Let us start by 
writing the consistency condition in the rest frame 

r- 1 <(1 + 2/ ) FdV) + AJ2 (pW) = . (54) 

If p = 0, then 
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H 3 > 



A oo + g A fe 



(l)* 



-r- 1 <(l + 2/ )FCi^)^0 



(55) 



and if p 



H) AS } = . (56) 
We see that A$ may be transverse but not both transverse and traceless. 

The point is that any term in A^J which is not strictly required by the consistency con- 
ditions may be absorbed into Affl, and so there is no loss of generality if we simply take 



with 



A 



A (i)ii = _ r -iA^ 
_ ((1 + 2/ ) Ftfwp) 



(57) 
(58) 



The second order equation then reads 



Q [FQ [C 2 ]] = -u 2 p A + (1 + 2/ ) F(%_ 



(59) 



Observe that the right hand side vanishes if integrated against u p , so the equation may be 
solved, but not if integrated against u> p / F. Thus the solution is 



C 2 = (l + 2 /o ) Cl 2 -A^-^ 



(60) 



where 



F 



<(i + 2 /o kc 2 }-a(^ 



(61) 



So far, we have 



Dpp p^p v 



|(l+2/ )C 1 2 -A^-^ 



(62) 



which is traceless. However, this cannot be the true correction to the energy-momentum 
tensor because H"!] ^ 0. This means that the parameter T in our equations is not the 
physical temperature T p h ys , but rather 
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T Tpfoy S ST 

where 5T is chosen so that 



d P, 5T = 4 P ^ = n°° 



dT T 

^ t- f~\ "v~\ n T * I ... 

2,phys 2,phys 



The physical correction to the energy-momentum tensor IT^ h = fl!] 4 h = and 



LL 2,phys — ii 2 g ii 2 

which is still traceless. 

Now there is only one traceless tensor quadratic in and so we must have 



where 



^ = 12r 2 ((l + 2 /o )^) 
It is interesting to recall the identities 

<(l + 2 /o )^)=T 2 A <u; 5 )=7T<a; 5 ) 

<(1 + 2/o) ^) = T 2 ± ( W 3) = 5T ( W 3) = 5T 3^ p = 15T 2 {p + p) 

Therefore for the Marie collision term F = T we get 



and for the Anderson and Witting collision term F = u p 25|, |26] 

K aw = 180t 2 T 2 (p + p) 
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In section [IV] we will see that this is precisely the form of obtained in the quadratic 
DTT. 

For clarity, we will briefly summarize the main logical steps followed in this section. 
For a linear collision term the variational equation ( l24"j) becomes (1341) . If we expand the 
Lagrange multiplier A MI/ and the nonequilibrium correction to the distribution function (now 
parametrized by the new variable ( given by ([2"01 to satisfy the H theorem), we can solve 
the variational equation perturbately up to second order in the relaxation time. The first 
and second order solutions are given by (I4"9j) and (1601 . respectively. We emphasize that the 
assuptions made regarding the integrability conditions, namely Aqq = Xf = Aq° = and 
that leading to (|57|) . constitute no loss of generality in the development. 



III. COMPARISON TO GRAD'S ANSATZ 



We will now compare, in the context of the method of moments, the closure provided by 
ihe entropy production variational equation with that provided by Grad's quadratic ansatz 
13] (see also 10. 11|; for a recent generalization to multicomponent systems see js]). 



Taking moments of the Boltzmann equation one obtains an infinite set of equations: 

d, J Dpp^f = = J Dp I col [f] = 
d,J Dpp»p»f = T^ = J Dpp»I col [f)=0 
dj Dpp^p u p s f = j Dpp u p s I col [f] 



(72) 



where in the first and second lines we have used that are collisional invariants. This 

infinite set is completely equivalent to Boltzmann equation. 

The method of moments rests on the assumption that a finite subset of this hierarchy 
will give a reasonable description of the hydrodynamic regime. The most common case is to 
consider only the first three equations of the hierarchy. However, the truncated system is not 
a closed one, since the derivative of the second moment can not be expressed solely in terms 
of the hydrodynamic variables (N fl ,T flu ). In order to close the system, and following Grad's 

14 
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idea, Israel and Stewart |9| proposed expanding the single-particle distribution / = / + 5f 
around its equilibrium value in a Taylor-like series in p M and truncating it at quadratic order. 
In the notation of eq. (jHJ) above, this is 

X = C+C^ + C^p„ + ... (73) 

where the coefficients (C, C M ) correspond to a shift in chemical potential, temperature and 
velocity, and may be taken as zero. The coefficient C^ u is subject to the constraint that we 
must reproduce eq. It therefore has to be traceless and transverse. In the rest frame, 
the nontrivial components satisfy 

2 (C ij co*) = U ij . (74) 

This constraint does not determine C % i by itself, and so we must resort to eq. (1721) or else to 
attempt a solution to the Boltzmann equation. The first two equations in (172]) . representing 
particle number and energy-momentum conservation, must be supplemented by an evolution 
equation for the dissipative tensor IP". The traditional way of obtaining the evolution 
equation is to use (173]) in the third equation of (1721 . In this way the Israel-Stewart equations 
are obtained. Recently, Denicol et al Q] reobtained these equations directly from the kinetic 
theory definition of the derivative of IP 1 ' instead of relying on the second moment equation. 
We note that the form of the equations obtained by these authors is the same as that of 
Israel-Stewart, but with different transport coefficients which result in better agreement with 

Boltzmann equation. In any case, a detailed analysis of the transport equation is necessary. 

r~i 

Luzum and Ollitrault j4| point out that the majority of works in this area assume that C % i 
is p-independent, therefore C 3 ' = IP/2 (u)p)- Comparingto eq. (|4"9]1 above, this corresponds 
to the case where F = T. A more detailed analysis (14, 5J) shows that while this obtains in 
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some cases, such as a A0 4 theory ([2|, |5j), it is not a good description of a hot gluon plasma. 
Moreover, the usual analysis also assumes for IP a gradient expansion 



IP = -ri<j ij + ... 
where a 1 ^ is the shear tensor, defined in the rest frame as 

15 



(75) 



a 



u' J + u 



3,i 



(76) 



and r) oc T 3 is the shear viscosity. 



On the other hand, 



also analyze more general options for C % K Concretely, they 



analyze cases where C* J oc uJ p a , with a = 0, 1/2 and 1. In the EPVM framework they 



correspond to F = T (co p /T) a . So a = corresponds to Mar 
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e's kinetic equation, while 



26|). 



a = 1 gives the equation proposed by Anderson and Witting 

In summary, we have found a way to associate a simple kinetic equation (which is nev- 
ertheless consistent with energy-momentum conservation and the Second Law) to the non- 
standard nonequilibrium corrections investigated in 4], 15]- Moreover, our treatment nowhere 
assumes a gradient expansion such as ( 175|) . We will discuss this point further in next Section. 

We believe that the form of \ we have found here (even at first order), being more flexible 
than the way Grad's ansatz is usually implemented, may be useful to improve the description 
of freeze-out in heavy-ion collisions, at least in a phenomenological fashion. 



IV. EPVM AS A DIVERGENCE-TYPE THEORY 



In this section we give a brief overview of divergence-type theories 
explicit the connection between DTTs and the EPVM. 



21. 
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A. Divergence type theories 

According to Geroch and Lindblom {2]]], the hydrodynamical description of a nonequi- 
librium state requires, besides the particle current and the stress-energy tensor T^, a 
new third order current A^ up . These currents are obtained as derivatives of a generating 
current x M with respect to the hydrodynamic variables a = fi/T, (3^ = u^/T and where 
fi is the chemical potential (which vanishes identically for a conformal theory), T is the tem- 
perature (see below), w M the velocity and the nonequilibrium tensor ^ v is symmetric and 
traceless. For a conformal theory ^ u is also transverse £^ v u v = 0. The relevant equations 
for a conformal theory are then 
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T» v = %- and (77) 



A^p = ^- . (78) 



The entropy flux is 



S» = x » - p v T^ - Cu P A^ p . (79) 
Assuming that T^ v is conserved, the entropy production is 

S% = -tvA™ (80) 

so knowledge of the entropy production gives an equation for A^ p and thus determines the 
evolution. 

Because T^ v is symmetric, x^ must derive from a potential 

X OP, { ° ' 

The most general conformally invariant potential containing up to quadratic terms in the 
nonequilibrium tensor and yielding a traceless energy-momentum tensor is 

X = aT 2 + T~% x u T u x - cT' 6 [&. A f * + 24£ T ^u V + 168 (f rA *x V) 2 ] . (82) 

This form of the potential assumes that has dimensions of T 4 . We retain terms involving 
£^ v u v , which are zero "on shell" but contribute to the derivatives of the potential. The 
derivatives are computed according to the rules 

cfT 

= T apn = tV , (83) 



dP v 



= TbT and (84) 
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whereby we get 



(85) 



X 11 = ^ 



|2aT 3 - 2T- 1 £ tX u t u x + 6cT~ 5 £ rA £ rA + 2A^{u T u x + 168 (£ T \u T u x )' 
+ 2A^ A {T-% x u r - 2AcT- 5 [^ e x u T + 14 (&*uV) fan*] } 



and, assuming a transverse £^ 



where 



and 



u u u^ + - A"" 
3 



+ 2£ 



48cT' 



4 tf tf rA 



p = 6aT 4 - 30cT~ 4 £ tA £ 



= T' 1 yg^ x u T + g^ T u x - -# T V 
- 24cT- 5 + C T u x - ^ tA w m J 



The entropy flux becomes 



i7) 



(89) 



5^ = 0" {8aT 4 - 36cT- 4 e rA r A } 



(90) 



We see from eq. ( 1881) that T is not the temperature as measured by a comoving observer. 
We have already encountered this situation in subsection III C 21 As in there the solution 
lies in a temperature shift (cfr. eq. fl63|) ). This redefinition of the temperature represents 
a correction of order £ 4 in eq. (157)1 and of order £ 3 in A^ uS , which we are going to neglect 
since they correspond to terms that would be obtained from a cubic generating function (see 



|20|). 
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B. To DTTs from EPVMs 



We will now show that the EPVM leads to the DTT discussed before. The idea is to seek 
a solution for ^ v as an expansion in the small parameter r. Therefore we write (compare 
toeq. flM}) 



+ t(2) 



(91) 



Matching the first order correction to T^ v in both theories we get 



tO)ij = 1_ (fp \ \(0)ij 
s 15 \ F 



(92) 



We use this result to compute the first nontrivial correction to the entropy flux and match 
to eq. (j52p . This determines the c coefficient 



5T 5 /u. 



12 



v 
F 



'<4 

F 2 



(93) 



Knowing c, we can math the full T^ v to get 



(2)ij 



K + —Tt 2 ( -f 
45 \F 2 



A (0), A (0 )mj _ 1^ A (0) A (0)mn 



(94) 



where F is defined in eq. (|6j 

The only remaining step is to find the equation of motion for ^ v by matching the entropy 
production as given in the DTT to the corresponding result from the EPVM, eq. ( |53l) . This 
gives (in the rest frame) 



C j = -r 



12 



r T 1 k 

5 T + 3^ 



ki 



where 



(95) 



■q = — — ( —) and 
' 15T \ F 



6 /^V" 1 /^ 



F 2 



(96) 
(97) 
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We recover eq. (I75|) when r — > 0. 



V. SUMMARY 

Relying on entropy production variational methods and using a linear collision operator 
that satisfies the H theorem, we have shown how to associate a one particle distribution 
function to the energy momentum and entropy currents of a conformal fluid, in a way that 
generalizes Grad's ansatz. The entropy production variational method leads to a nonequi- 
librium correction to the distribution function which at first order, and for a specific form 
of the linear collision operator, reproduces Grad's quadratic ansatz. For other choices of the 
collision operator we obtain a generalization of Grad's ansatz, in which the nonequilibrium 
distribution function can have a dependence on momentum other than quadratic as indicated 
by recent developments 0,0]. 

Moreover, by equating the entropy production to its extremal value we obtain a hydro- 
dynamics which belongs to the class of divergence-type theories. The kinetic origin of the 
DTT ensures that the Second Law is satisfied. The derivation provides a kinetic interpreta- 
tion of the DTT nonequilibrium tensor as a Lagrange multiplier enforcing the stress-energy 
constraint when extremizing the entropy production. 

Although based on a linear collision term, we believe that our results provide an interesting 
link between kinetic theory and dissipative divergence-type theories which may prove useful 
in the study of heavy-ion collisions, for instance, to improve the description of the freeze-out 
process. It would be interesting to find the dependence of the relaxation time on energy 
by matching the results of hydrodynamic simulations of heavy-ion collisions to data, in this 
way constraining the form of the linear collision operator used in the kinetic description of 
relativistic plasmas. 
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